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Abstract 

This work examines the large-scale deployment of energy harvesting sensors for the purpose of 
sensing and reconstruction of a spatially correlated Gaussian random field. The sensors are powered 
solely by energy harvested from the environment and are deployed randomly according to a spatially non- 
homogeneous Poisson point process whose density depends on the energy arrival statistics at different 
locations. Random deployment is suitable for applications that require deployment over a wide and/or 
hostile area. During an observation period, each sensor takes a local sample of the random field and 
reports the data to the closest data-gathering node if sufficient energy is available for transmission. 
The realization of the random field is then reconstructed at the fusion center based on the reported 
sensor measurements. For the purpose of field reconstruction, the sensors should, on the one hand, be 
more spread out over the field to gather more informative samples, but should, on the other hand, be 
more concentrated at locations with high energy arrival rates or large channel gains toward the closest 
data-gathering node. This tradeoff is exploited in the optimization of the random sensor deployment 
in both analog and digital forwarding systems. More specifically, given the statistics of the energy 
arrival at different locations and a constraint on the average number of sensors, the spatially-dependent 
sensor density and the energy-aware transmission policy at the sensors are determined for both cases 
by minimizing an upper bound on the average mean-square reconstruction error. The efficacy of the 
proposed schemes are demonstrated through numerical simulations. 


L Introduction 

Wireless sensor networks (WSNs) eonsist of spatially distributed sensors that have the ability 
to sense the physieal environment, proeess the gathered information, and eommunieate through 
the wireless interface. In recent years, WSNs have been adopted in a wide range of applications, 
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such as environmental monitoring, disaster reeovery, and battlefield surveillance, ete [[U. In 
these applieations, sensors are often deployed in large-seale and in hostile areas making human 
maintenanee and battery-replaeement impraetieal. Due to these reasons, energy harvesting teeh- 
niques for sensor nodes O-JH, which enable the conversion of ambient energy (such as solar 
[[5l , vibrational Q, or thermal energy |I71) to eleetrie energy, are being developed and used to 
prolong sensor lifetime. By employing energy harvesting teehnology, the characteristics of the 
energy arrival and the effieieney of energy usage will have a signifieant impaet on the sensing 
performanee. It is therefore neeessary to adapt the sensor deployment and sensor operations to 
spatial variations of the energy arrival proeess. 

The main objeetive of this work is to determine optimal sensor deployment strategies, namely, 
spatial densities of energy harvesting sensors, for the purpose of sensing and reeonstruction of 
a spatially eorrelated Gaussian random field. The sensors are assumed to be deployed randomly 
and in large scale aeeording to a spatially non-homogeneous Poisson point proeess (NHPPP) 
if^- lfTOl . During an observation period, eaeh sensor takes a loeal sample of the random field and 
reports the observation to the closest data-gathering node based on a threshold-based energy- 
aware transmission eontrol policy. The policy allows the sensor to transmit only if it has 
aeeumulated enough energy for transmission. We eonsider both analog-forwarding (AF) systems, 
where sensors transmit a scaled version of their analog measurements to the data-gathering node, 
and digital-forwarding (DF) systems, where digital representations of their measurements are 
forwarded instead. The random field is then reeonstrueted at the fusion eenter based on the 
information gathered from the sensors. 

In this work, we assume that the sensors’ operations are supported solely by ambient energy 
(i.e., energy harvested from the environment) and, thus, their transmit powers and probabilities 
depend strongly on the characteristics of their energy arrival at their respective loeations. To 
reduee the field reeonstruction error, the sensors should, on the one hand, be more spread out 
over the field to gather more informative samples, but should, on the other hand, be eoneentrated 
more at loeations with large energy arrivals or with large ehannel gains toward the closest data- 
gathering node. This tradeoff is exploited to determine the optimal random sensor deployment in 
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both AF and DF systems. In particular, given the locations of data-gathering node^li the statistics 
of the energy arrival at different locations and a constraint on the average number of sensors, we 
determine the spatially-dependent sensor densities and the energy-aware transmission thresholds 
at the sensors by minimizing an upper bound on the average mean-square reconstruction error. 
The energy-aware transmission policy allows a sensor to transmit only when its accumulated 
energy (i.e., its battery level) is beyond a certain threshold and remains silent, otherwise, so 
that energy can be preserved for use in later time slots. Notice that, in this work, we consider 
only a simple threshold-based transmission policy and focus on the global effect of the sensor 
deployment problem. In general, the optimal energy-aware transmission policy may involve 
continuous power control in accordance with the data traffic, the transmission deadline, and queue 
stability etc, but is not considered in this work to maintain tractability. Readers are referred to 
ifm - lfTBll for further studies on this topic. The efficacy of the proposed schemes are demonstrated 
by numerical simulations. 

In the past, sensor deployment problems have been examined mostly for sensors without 
energy harvesting capabilities (see [[T4ll for a survey on these topics). In particular, without 
energy harvesting considerations, sensor deployment policies have been proposed with the goal 
of minimizing the field reconstruction error in ifTSll . of guaranteeing connectivity in [fT^ . and 
of maximizing sensor coverage in ifTTll . In these works, sensors were assumed to be placed 
at deterministic locations, in which case, the task of finding the optimal sensor placement is 
often NP hard lfT4]| . Therefore, heuristic or approximate solutions were proposed to reduce 
the computational complexity. Due to advances in energy harvesting technology (H-im, similar 
problems have also been examined recently, e.g., in ifT^ and ifT^ . for energy harvesting wireless 
systems with considerations on the stochastic nature of the energy arrival at each node. In 
particular, in ifT^ and ifT^ . the deployment of energy harvesting relay nodes in sensor networks 
were examined with the goal of enhancing network throughput and of guaranteeing connectivity, 
respectively. However, they did not consider the cross-layer impact of the spatial dependencies 
of the energy arrivals and sensor measurements on the field reconstruction performance. In 

'Here, we assume that the locations of the data-gathering nodes are fixed and focus on the deployment of energy harvesting 
sensors. This is often the case in practice since data-gathering nodes are typically grid-connected and, thus, their placement can 
he more restrictive. However, in certain cases, further flexibility may also be given to the deployment of data-gathering nodes. 
Readers are referred to fini for further discussions on this topic. 
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our work, we consider the optimal random sensor deployment strategy for field reconstruction 
by taking into consideration the spatial correlation between the sensor measurements and the 
energy arrivals. Rather than designing a deterministic deployment scheme which places sensors 
at precise locations, we consider the random deployment of sensors and determine the optimal 
sensor densities at different locations. We argue that random deployment ll2T]| is more 
practical when the number of sensors to be deployed is large or when the sensors are to be 
deployed in a hostile environment. Comparisons between deterministic and random deployment 
schemes can also be found in ll22ll and [|2^ . 

Field reconstruction and decentralized parameter estimation are essential applications of WSNs 
and have been studied extensively in the literature for WSNs without energy harvesting capabil¬ 
ities, e.g., in [I24l - [|2^ . For given sensor locations, these works focused on the design of sensor 
transmission schemes under different centralized fusion rules. Two transmission systems have 
been considered the most in the literature, namely, AF If2^ - lf2^ and DF lfT7l - ll29l systems. In 
these systems, the amplifying gains [|2^ - [I^ and the number of quantization bits [I27ll - ll29l used 
for transmission by the sensors respectively can be chosen to minimize the field reconstruction 
error. Both of AF and DF systems are examined in our work under additional energy harvesting 
considerations. Decentralized estimation under energy harvesting constraints have also been 
investigated recently in [[^ and ifSTI . However, the sensor deployment problem and the impact of 
the location-dependent sensor measurements and energy arrivals on the estimation performance 
have not been explored before. A preliminary version of our work can be found in ll^ . but 
only for the AF case. 

The rest of this paper is organized as follows. A general description of the system model is 
presented in Section |nl The sensor deployment and transmission control policies are then derived 
separately for AF and DF systems in Sections |nl] and UVl respectively. Numerical results are 
provided in Section 0 and some concluding remarks are given in Section ??. 

II. System Model and Problem Deeinition 

Let us consider a WSN that consists of a large number of sensors deployed randomly in the 
region V according to a spatially NHPPP ifTOll with deterministic densities {A(v)}vev, where A(v) 
is the density at location v G V. In an observation period, each sensor takes a local sample of the 
random field and forwards it to the closest data-gathering node. The locations of data-gathering 
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Fig. 1. Illustration of the sensor deployment, the energy arrival distribution and the field values for a specific t. 


nodes can be arbitrary but are assumed to be fixed over time. The random field is denoted by 
{a:(v)[t]}vev> where x(v)[t] is the value of the field at location v G V in the t-th observation 
period. Following ifTSll and ^2M . we assume that the set of random field values {x(v)[f]}vev 
is jointly Gaussian with zero mean and with covariances E[x{\'i)[t]x{vj)[t]] = 
for all Vj,Vj G V, where al is the variance of x(v)[f], for all v G V, and p(vj,Vj)[f] is the 
correlation coefficient between x{vi)[t] and x{vj)[t]. The random field is also assumed to be 
stationary over time. Moreover, let {a(v)[f]}vev be the set of random energy arrivals, where 
a(v)[f] is the energy arrival at location v G V in the f-th observation period. An example of 
the sensor deployment, the field values, and the energy arrivals for a specific t are depicted in 
Fig. [U In practice, the set of sensor densities {A(v)}vev to be derived in this work provides a 
guideline for the numbers of sensors that should be scattered at different locations, e.g., from 
an air vehicle. 

To simplify our computations, let us partition the region of interest, i.e., V, into M disjoint 
subregions Vi,..., Vm with equal size. That is, we choose Vi,..., Vm such that V = V* 
with areas |Vi| = A, for all i, and with Vif]Vj = 4>, for all i ^ j. We assume that A is 
sufficiently small (compared to the variations of the random field) so that the sensor densities 
remain approximately constant in each subregion (i.e., A(v) Aj, for all v G Vi and for all i). 
Then, by assuming that AjA <C 1, the probability that exactly one sensor exists in subregion V* 
can be approximated as ^ ^ probability that more than one sensor exists in 

Vi is negligible. Moreover, with A sufficiently small, we can also approximate the field value 
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and the energy arrival at a sensor in subregion Vi by their respeetive values at the eenter of 
the subregion, whieh are denoted by Xi[t] and ajf], respeetively. Therefore, if a sensor exists in 
subregion Vi, its loeal observation ean be expressed as 


Xi[t] = Xi[t] + ni[t], 


( 1 ) 


where ni[t] is the observation noise at the sensor loeated at the subregion Vi. The observation 
noise is assumed to be independent and identieally distributed (i.i.d.) aeross sensors and over 
time aeeording to a Gaussian distribution with mean zero and varianee denoted by AA(0, cr^). 

To preserve energy at the sensors, we adopt a threshold-based energy-aware transmission 
poliey where a sensor transmits only if its aeeumulated energy exeeeds a eertain threshold, 
but exhausts all its energy when doing so. This seheme is referred to in the literature as the 
integrate-and-fire [I33]|, 041 transmission poliey. The energy threshold at the sensor in subregion 
Vi is denoted by 7 * and is ehosen to minimize the field reeonstruetion error in later seetions. 
The aeeumulated energy at eaeh sensor, say, the sensor in subregion Vi, varies over time t and 
ean be expressed as 



ei[t - 1] + ai[t], ifej[f-l]<7i 

ai[t], if ei[t - 1] > 7i- 


( 2 ) 


at the beginning of the f-th observation period. By assuming that the proeess {ej[f]}“o 
stationary over time, the probability that the sensor transmits in a given time slot is F^.^ji) = 
Pr(ei[f] > 7 j), whieh is the eomplementary eumulative density funetion (eedf) of ei[t]. An 
example of the ease where {ei[f]}“o stationary is given in the following seetions by eonsidering 
a Bernoulli energy arrival proeess. By the stationarity of the random field {xi[t]}'^Q and of the 
aeeumulated energy proeess {ei[t]}“o’ P we shall omit the time index t in later diseussions 
and foeus on the field reeonstruetion performanee at a partieular instant in time. 

Let us define a random variable Oi to indieate the presenee or absenee of a signal transmitted 
by a sensor in subregion Vi. Speeifieally, we set Oj = 1, if a sensor exists in subregion Vi 
and has energy above the threshold 7 ^, and set o* = 0, otherwise. In this ease, Oj ean be 
viewed as a Bernoulli random variable with probability Pr(oi = 1) = XiAFe^{'Ji) and Pr(oj = 
0) = 1 — AiAFei( 7 i)- Moreover, let rj be the signal reeeived from the sensor in subregion Vi 
at the elosest data-gathering node over a noisy ehannel. When Oj = 1, the reeeived signal rj 
depends on the the sensor’s observation Xi, the aeeumulated energy Cj, the type of transmission 
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scheme, and the quality of the ehannel towards the elosest data-gathering node. When o, = 0, 
nothing is reeeived and, thus, we set = null. Two types of sensor transmission sehemes are 
eonsidered, namely, the AF and DF sehemes. In the AF seheme, eaeh sensor simply transmits 
an amplified version of its reeeived signal to the elosest data-gathering node. In the DF seheme, 
eaeh sensor first quantizes its measurement into a binary representation veetor and then forwards 
it to the elosest data-gathering node, where the binary veetor is deeoded and the quantized sensor 
measurement is reeonstrueted. 

The signals reeeived at the data-gathering nodes are transmitted over the baekhaul to the fusion 
eenter where an estimate of the field values {a:(v)}vev are eomputed. Let o = [oi,... ,om]^ 
be the veetor representing the sensors’ transmission status. The transmissions from the data- 
gathering nodes to the fusion eenter are assumed to be error-free and, thus, are also the 

observations available at the fusion eenter. By assuming that the fusion eenter is aware of the 
sensors’ loeations and transmission statuJl i.e., it extraets only the set of field-bearing 

observations {ri\l < i < M, o* = 1} for field reconstruetion. The linear minimum mean square 
error (LMMSE) estimator is then adopted to obtain an estimate of eaeh point in the field using 
these observations. Interestingly, we show in Appendix |A] that the LMMSE estimate obtained 
with the set of observations {ri|l < i < M,Oi = 1} is equivalent to that obtained with the 
effeetive reeeived signal veetor y = [ 7 / 1 ,..., where 

{ Vi, if Oi = 1 

0, if Oi = 0, 

for all i. Therefore, the estimate of the field value a;(v) at loeation v ean be written as 


i;(v) 


po po 
'^a:(v)y'^yy Ti 


(3) 


where C“, A B|yy^|o|, C;,.,, = E'[x(v)y^|o], and ^ represents the Moore-Penrose pseudo 
inverse. The resulting MSE of the estimate on x(v) is then given by 

5(v|o) = <7^ - tr (4) 


^Here, we assume that the fusion center is aware of the locations of the sensors once they have been deployed and the 
transmission status of the existing sensors. The former can be obtained through positioning techniques after deployment and 
the latter can be obtained by appending the sensors’ IDs to their transmissions or by performing signal detection at the fusion 
center. However, the proposed methodology can also be used for the case where no knowledge is assumed at the fusion center. 
In this case, the received signal in the absence of a sensor transmission in subregion Vi should be n = Wi. 


( 5 ) 


and the average MSE over the entire sensor field is defined as 

5^^(v|o)Pr(o)dv, 

where O is the set of all possible realizations of o. Notiee that the above expression involves the 
summation over all possible realizations of o, whieh ean be intraetable in praetiee. To simplify 
our eomputations, we eonsider instead an upper bound of the average MSE given as follows: 



[ K 

'vev 

min E — k^y P ol 

keRM LI V ; J1 1 j 

dw 

(6) 

1 

M 

/ min Eo \E I" a:(v) — k^y ^ oil dv 

L LI V ; 1 1 JJ 

(7) 

1 

M 

/ min E T x(v) — k'^ypl dv 

ivevkeR^ ^ ^ 


(8) 

1 

M 

/ 

tr dv A 

■Cupper 5 

(9) 


where Cx(v)y — E[x{v)y'^] and Cyy = E[yy'^] are eomputed by taking the expectation over all 
possible realizations of o, and thus, do not depend on the actual subset of transmitting sensors. 
The vector k inside the expectation in db]) can be viewed as the linear estimator that should 
be optimized separately for each given value of v and o. The inequality in (|7]) follows since 
the minimization over each term inside the expectation must be smaller than the minimization 
over the entire expectation (which amounts to using the same k for each o). The latter can be 
interpreted as the minimum MSE attainable when the fusion center does not have knowledge of 
the sensors’ locations and transmission status. This bound is tight only when the optimal linear 
estimator k is approximately the same for all realizations of o that occur with high probability. 
Even though the sensors are deployed randomly, their locations are fixed once they are deployed 
and do not change over time. The MSE in (jb]) provides a measure of the average performance 
over all possible sets of sensor locations following an NHPPP with sensor densities but 

is only an approximation of the actual MSE in practice, which corresponds to just one realization 
of o. 

The main objective of this work is to determine the optimal sensor densities and 

the energy thresholds { 7 *}^! by minimizing the MSE upper bound in ([91) for both the AE and 
the DE systems. The optimized sensor densities should achieve the optimal balance between the 
energy arrival probability, the channel gain, and the sensor field correlation whereas the optimized 





Fig. 2. Illustration of filed reconstruction for the analog forwarding system. 


energy thresholds should exploit the tradeoff between the sensor’s transmission probability and 
the reception quality at the data-gathering nodes. The two systems are examined separately in 
the following sections. 


III. Optimized Sensor Densities and Energy Thresholds for 
Analog-Forwarding Systems 


In this section, we determine the optimal sensor densities and energy thresholds 

in AF systems based on the minimization of the average MSE upper bound in @. 

In the AF system, each sensor transmits an amplified version of its local observation to the 
closest data-gathering node, as illustrated in Fig. [21 In this case, the signal transmitted by a 
sensor in subregion Vi can be written as 


•y/NjCj 


s, = ^=^Xi, ( 10 ) 

where n* > 0 is the amplifying factor, is the energy available at the sensor in Vi, and the 
received signal at the closest data-gathering node is 

\/ Nje* 


yi= [hi 


■.Xi+Wi 


V^x + 

where hi is the channel coefficient between the sensor and its closest data-gathering node 
and Wi ~ AA(0, is the additive white Gaussian noise (AWGN). Here, we do not consider 


( 11 ) 

i<fl 


^The signal model under consideration can also accommodate broadcast transmissions, where the signal transmitted by a 
sensor in subregion Vi may be received simultaneously by multiple data-gathering nodes simultaneously. In this case, hi can he 
viewed as the effective channel coefficient after signal combining at the fusion center. 
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the effects of short-term fading and assume that hi is a constant that depends only on the 
distance between the sensor and its data-gathering node. Therefore, by defining gi = 
as the effective channel gain experienced by the sensor in Vi, the received signal is expressed 
as yi = {giXi + Wi)oi. For simplicity, we assume that the sensor depletes its battery every time 
it transmits (i.e., any residual energy is dumped after each transmission). Notice that, since the 
observation Xj is random, the energy required to transmit the signal in (fTOl) may exceed the 
accumulated energy e* with a certain probability. In this case, a saturation effect may occur, 
causing additional distortion on the received signal. However, by choosing Ki to be sufficiently 
small, the probability that the saturation occurs is small and, thus, is omitted for simplicity. 
Studies on the impact of this effect on the distributed estimation performance is beyond the 
scope of this work, but can be found in ll^ . 

For convenience, let = Pr(oj = 1) = Fe^{'yi}XiA be the probability that a sensor exists 
in Vi and transmits in the given observation period. Moreover, let gii'ji) = E[gi{'yi)\ei > 7 ^] 
and gfi'ji) = E[gf{'yi)\ei > 7 *] be the conditional first and second moments of the effective 
channel gain given that the sensor transmits. In this case, we have Ca;(v)y = C 3 ;(v)xDgDa> where 
Dg = diag (^ 1 ( 71 ),..., (?m( 7 m)) and Dq, = diag(ai,..., om); thus, the MSB upper bound in 
dH) can be written as 

Supper,AF = - tr ($DgD„CyyD„Dg) , (12) 


where ^ - p /vgv Cx(v)xCx(v)xC?v is defined such that 

(T^ f 

Wi,j = HT / Vi)p(v, Vj)dv = 

I Jv£V 

for all i,j, and the (*, j)-th element of Cyy can be derived as 


(13) 




(14) 


for i=j 
for 

for all i,j. Notice that $ is a matrix that depends only on the correlation of the sensor field and 
not on the optimizing parameters {Xi}fii and { 7 *}^^. Let us take the eigenvalue decomposition 
of $ so that $ = USU^, where S is a diagonal matrix consisting of the eigenvalues of $ and 
U is a unitary matrix consisting of the corresponding eigenvectors. In this case, the term in (fT^ 
can be written as 


e 


upper 


,af - cr. 


- tr (^S5U^DflD„C-!D„D,US^ 


'g^a^yy^a-^g 


(15) 
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To further facilitate our derivations, we introduce the following lemma from [|25]| . which is a 
consequence of the Cauchy-Schwarz inequality. 


Lemma 1 ( llISll l. For any G E and positive-definite matrix Q G the following 

inequality holds: 


tr(G^Q-iG) > 


[tr(G^G)]' 
tr(G^QG) ■ 


(16) 


The equality holds when G^ Q 2 = cG^ Q2, where c is a constant. 


By taking G = DQ,DgUS 2 and Q = Cyy, it follows that 


e 


upper,AF ^ ^ 


tr ( StU'^YDgD«) 2 USt 


tr (SsU^D^D^CyyD^DgUSs 


2 m 2 


[tr($(DgDj^)] 

^ tr($D3D„CyyD„D§) 


'bobj,AF('^5 T); 


(17) 

(18) 


where A = [Ai,...,Am]^ is the vector of sensor densities and 7 = [ 71 ,..., 7 ^]^ is the 
vector of energy thresholds. By Lemma [U the bound is tight when (DcD^US = 

c(DQ,DgUS2)'^Cyy^, for some constant c. The upper bound ^obj,AF(A, 7 ) is then utilized as the 
objective function for optimizing A and 7 in the AF case. Even though the bound may not 
be tight in general, it captures the essential behaviors of the MSE with respect to the sensor 
densities and transmission thresholds and, thus, allow us to obtain a tractable solution to an 
otherwise intractable problem. The effectiveness of the solution is demonstrated in Section |Vl 
Specifically, let us consider the optimization problem where the MSE upper bound .^obj.AF in 
(fTSl) is minimized subject to a constraint on the average total number of sensors, i.e.. 


min 

^,7 

'Cobj,AF(A, 7) 

(19a) 


M 


subject to 

5^A,A<A, 

i=l 

(19b) 


0 < AjA< ca, for z = 1,..., M. 

(19c) 


where A is the constraint on the average number of sensors in the network. Here, the approximate 
probability that a sensor exists in subregion Vi, i.e., AjA, is restricted within [0, ef\ and ca < 1 
should be small enough to make the approximation accurate enough. The value of A* is assumed 







12 


to be positive to ensure that Cyy is invertible. Following the relation in Q, the energy threshold 
7 determines the distribution of the accumulated energy process and, thus, affects the objective 
function through the probabilities a*, for all i, which involve the event that a sensor exists and 
transmits in each subregion. 

Since the first term in (fTSl) is a constant, the optimization problem can be equivalently 
formulated as the maximization over the second term in (fTSl) . By the change of variable A* = A* A 
and by defining A = [Ai,..., A^]^, the optimization problem can be written explicitly as 
follows: 


max 

A,7 

Jaf(A,7) = 

4f(A,7) 

(20a) 


M 



subject to 

5^A,< A 

2 = 1 


(20b) 


0 < Ai< ca, 

for f = 1 ,..., M. 

(20c) 


The term in the numerator of the objective function is 


M 




. i=l 


and that in the denominator is 


M 


2=1 

M 


4f(A,7)=5^0M + 




( 21 ) 


( 22 ) 


where a? = 

Notice that the solution to the optimization problem in (l20l) depends on the statistics of the 
energy arrival. In the following, we provide an example of how such a problem can be solved 
under the Bernoulli energy arrival model. 


A. Solution for Bernoulli Energy Arrival Case 

In this subsection, we apply the above techniques to cases where the energy arrival can be 
modeled as an i.i.d. Bernoulli process. This model has been widely adopted in the literature, 
e.g., in [[3^ - [[39l . due to its tractability and because of its ability to model sources with sporadic 
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1 - Pi 



Fig. 3. Energy harvesting model of Bernoulli energy arrival. 


energy arrivals in discrete-time. Examples of such energy sources may include vibrational, RF 
energy, or wind etc. It can also be used to approximate more stable energy sources, such as solar, 
by choosing each time slot to be sufficiently small so that energy will arrive approximately at a 
constant rate when viewed from a larger time-scale. Other energy arrival models can also fit into 
our framework, but may require different optimization techniques from that to be introduced in 
the following. 

Suppose that the energy arrival process consists of a sequence of i.i.d. Bernoulli 

random variables with Pr(aj[t] = 5) = pi and Pr(aj[t] = 0) = 1 — pi, Vf. That is, the sensor in 
Vi harvests energy 5 in each observation period with probability pi and fails to harvest energy 
with probability 1 — Pi. Since the sensor expends all its energy once the energy threshold 7j is 
exceedeqd, the accumulated energy forms a (iVj-|-l) -state Markov process, as illustrated 

in Fig. [3l where iV* = min{?7, : n6 > 7*}. The stationary (or steady-state) distribution is given 
by Pr(ei[f] = 0) = Pr(ei[t] = Ni6) = and Pr(ei[f] = /c(5) = ^ for fc = 1 ,..., - 1 . 
In this case, the probability that the sensor in subregion V* transmits (if it exists) is given by 
Feiili) = Pr(ei[f] = Ni5). Notice that, in the Bernoulli energy arrival model (c.f. Fig. [3]), the 
accumulated energies are integer multiples of 5 and, thus, the energy threshold 7^ can also be 
set as a multiple of 5, in which case, we have 7* = Ni5. Consequently, we have Fe.{'yi) = 

and where hi = y/nihi. Notice that the stationary distribution as 

well as these values change as the energy threshold 7* is adjusted. By substituting the above 
into (1^ and (l22l) . we get 


M 


Tnum 

^AF 


(A, 7) ={Y1 

i=l 


(23) 


''The accumulated energy may not be depleted after each transmission in the AF case since the required transmit power is 
random. However, we assume for simplicity that the remaining energy is omitted after each transmission. 
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and 


M 




2=1 


M 


i=l j^i 

and, thus, the optimization problem in (l20l) can be written as 

M M \2\2 


h~'- 

] h 


max —— 

A ,7 M 


i=l j=l 


lilj 


A3 ^ A3 ^ 


,, ,, 

ECiA+E^'i+EE^y— 

i=l j^i 


tr tEi '^77, 


M 


subject to ^ Aj < A, 


i=l 


0 < Ai< ca, i = 1,...,M, 

7 j = Zj( 5 , Zj G N, i = 1 ,..., M, 


(24) 


(25a) 


(25b) 

(25c) 

(25d) 


where Aij = (j)i^i(j)jjh^hjp^p^jd'^ , C* = (pi^JijpfS^al , A = (pi^ia^p^^d'^hfal , and (7*,^- = 

A,j(ylpiyi, ^j)pipfhTh]5^- 

Notice that this problem is a mixed integer nonlinear programming (MINLP) problem, which 
is difficult to solve in general. To address this issue, we consider a relaxation where the integer 
constraint on 7* is replaced with the inequality constraint 7* > 5. Then, by defining f{A,j) = 
J2iLi the optimization problem can be written equivalently as 


^ AS ^ A3 ^ 




mm 

A ,7 


7i ^ Ti 111 3 

2=1 2=1 ^ 2 = 1 J^2 


/(A, 7) 


M 


subject to ^ Aj < A, 


2=1 


0 < Ai< ca and 'ji> 5 , for z = 1,..., M. 


(26a) 

(26b) 

(26c) 


Notice that the relaxed problem is still non-convex, but can be approximated by a series of 
geometric programming (GP) problems using the condensation method ||40| . 

Specifically, in each iteration of the condensation method, the function /(A, 7) is replaced 
by its monomial approximation so that the objective function in (I26al) becomes a posynomial 
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Algorithm 1 Optimized Sensor Densities and Optimized Energy Thresholds (OTOET) Seheme 
in the AE System 

Initialization: Set £ = 0 and a solution accuraey e > 0. Eind a feasible solution 

Iteration: 

(i) Compute (3^^ , for j = 1,..., M, by (|2^ . 

(ii) Replaee /(A, 7) with /*'^^(A, 7) and solve (l2^ . Eet the solution be j. 

(iii) Repeat (i) and (ii) until 

JAF(A''’,f')) 


and that the problem in (|2^ ean be formulated as a standard GP problem. More speeifieally, 
let be the solution obtained in the £-th iteration of the condensation method. Then, 

based on the inequality between arithmetic and geometric means iHTll . it follows that 

A.jA;A; \gS 

A/‘'>(A,7), (27) 


M M 


M M 


A. .A^AS 

/(A.7)=y:y:^“>nn 


i=l 7=1 


7*7i 


2=1 7 = 1 


where 


qW _ A , 
Pi,j ~ 


(Al'')=(A? 



/(a''>.7> 


(28) 


By approximating /(A, 7) with p^\A, 7), the optimization problem in (|2^ becomes a standard 
GP, which can be converted to a convex optimization problem and solved using the interior point 
method [|42l|. Note that solving the problem given in (|2^ by replacing /(A, 7) with 
yields a solution that is also feasible in (l26l) and can be used to find the solution in the next 
iteration of the condensation method. By initiating with a feasible solution (A^°\ 7^°^), the above 
process can be repeated until the objective value converges. The procedure is summarized in 
Algorithm [TJ An approximated solution to the original problem in (1251) is thus obtained and the 
optimized energy threshold is further rounded to the nearest multiple of 5. The effectiveness of 
this scheme is demonstrated through Monte Carlo simulations in Section |Vl 
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Following [|42l Chapter 11], the GP problem in (|2^ ean be eonverted to a eonvex optimization 
problem with n = + 4M parameters and m = 2M^ + 6M + 1 inequality eonstraints. The 

objeetive and the eonstraints of this problem ean be further used to synthesize a log-barrier 
funetion that satisfies the self-eoneordant property, and Newton’s method ean then be used to 
solve it. The eomplexity of eaeh Newton step grows as 0{mm? + n^) [|4^ and the number of 
Newton steps required ean be bounded by y/m ll42l . Henee, the eomputational eomplexity of 
the problem grows as Details are omitted due to spaee limitations, but ean be obtained 

following the steps in ll4Sll . 


B. Complexity Reduction 


It is neeessary to note that, even though the above approaeh ean yield good solutions to the 
random sensor deployment problem, the eomplexity ean be high when the number of subregions, 
M, is large. This is not a problem in most oases sinee the sensor densities and energy thresholds 
need only be oomputed offline. However, if it is neeessary to reduoe the eomplexity without 
redueing the resolution of the subregions, one ean further reduoe the number of parameters 
by assuming that the parameters in neighboring subregions are the same. This assumption is 
reasonable sinee the dimensions of a subregion is assumed to be muoh smaller than the spatial 


variations of the sensor field and the energy arrivals. 

Suppose that the M subregions are oombined into N olusters Ci, C2 ,..., Cat, eaeh oonsisting of 
Me = M/N subregions. Let Ac = [A(i), A(2 ),..., A(jv)] and 7^- = [7(i), 7(2), • • •, 7(JV)], where 
A(„) and 7(„) are the values of A* and 7*, respeotively, for all i G C„. Then, by letting A(^n),{m) — 

~ ^ A) ~ ^ A) ~ '^iGCn j &Cm.,i¥=j 


the optimization problem in S26\i ean be reduoed 


mm 

Mac 


N 

E 

n=l 


A3 

(”■) I n (”') 

e^(n) Tz^(n) 


7(n) 


/ G(n),(m) 

A(n)7(r.) 


\nY\m) 


A(n) 


/c(Ac,7c) 


(29a) 


N 

subjeot to ^ McA{n) < A, (29b) 

n=l 


0 < A(„)< eA, 7(n) > 5 , n = l,...,N. 


(29e) 


The problem ean then be solved using the oondensation method, similar to that in Algorithm [H 
and the eomplexity is redueed to 0((M/Mc)^). 
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Fig. 4. Illustration of filed reconstruction for the digital forwarding system. 


IV. Optimized Sensor Densities and Energy Thresholds for Digital 

Forwarding Systems 

In this section, we optimize the sensor densities the energy thresholds {'jijfii in 

DF systems based on the minimization of the average MSE upper bound in ®. 

In the DF system, each sensor first quantizes its local measurement into a binary representation 
vector and forwards it to the destination, where the measurement is reconstructed. An illustration 
of the DF system is given in Fig. HI Suppose that Bi is the number of quantization bits at a sensor 
in subregion Vj. The sensor measurements at the sensor are quantized into 2^* representation 
levels given by the set Aii = {rrii^i,... using the quantization function Qi defined by 

Qi{xi) = min{mi E M.i ; \xi — mi\ < \xi — m'|, Vm' E Aii}. The index of the representation 
level is then encoded into the binary vector bj and transmitted to the closest data-gathering 
node. Suppose that mi = Qi(xi) is the quantized value of Xi at the sensor in Vi and mi is the 
corresponding value reconstructed at the data-gathering node based on its received signal. In 
this case, the effective received signal at the data-gathering node can be written as yi = rhiOi 
and, together with the quantized measurements from other sensors, is utilized to perform the 
FMMSE estimate at the fusion center. The MSE and its upper bound can be written similarly 
as ([5]) and dH), respectively. Notice that the MSE is affected by both the quantization error that 
appears when representing Xi with mi and the channel error that causes the difference between 
mi and m^. 

Following the procedure in ll27]l . we adopt a uniform quantizer, where mi^i = ^ ^ ^ 
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for / = 1..., 2^% and Aq. = for W chosen sufficiently large such that Pr(|xj| > W) ~ 

0. In this case, the quantization error €i = rrii — xi is bounded in and can be 

approximated as a uniform random variable over this region when Aq. is sufficiently small. In 
this case, the variance of e, is given by Moreover, it has been shown in ETII and 

lf4^ that, for Bi sufficiently large, for all i, the quantization errors are approximately 

uncorrelated and independent of the sensor measurements Once the quantized value 

rrii is obtained, it is converted into a Bi-hit vector bj using a natural binary code (where b* 
is taken as the binary representation of the quantization index) and transmitted to the data- 
gathering node. By adopting BPSK modulation ||271 . the bit-error probability of the transmission 
is ei = Q{\ -3-^), where hi is the channel coefficient to the closest data-gathering node, Ci/Bi 

2 /V, .-S ^ 

is the energy per bit, and Q{u) = e~^du. By assuming that bj = ..., 4 is the 

binary vector received by the data-gathering node, the reconstructed quantization level is then 
rhi = [2(X]f=i -f 1) — 1 — since J2k=i -f 1 is the integer value of the 

binary vector bj. Properties of the quantization and channel errors are utilized to obtain explicit 
expressions of the MSB upper bound. 

In obtaining a tractable upper bound, let y A where yi = rriiOi, be the vector 

of received signals at the data-gathering node when the channel is noiseless. By also letting 
y = [^1 ,... iVivlY, where Vi ^ Vi — Vi, the MSB upper bound in ([8]) can be further upper- 
bounded as 


e 


upper,DF 


|V| 



min B'flxfv) — k^y -(- k'^y 

keiR^f 


k^y^dv 


2 f 

— / mm 


{E [|x(v) - k^y|2] + E [|k^y - k^yl^] } dw 


2 _ (tr (^DaDg))^ 

" tr($D«(Cyy + Cyy)D„) 


(30) 

(31) 

(32) 


A 


^objjDF 


(33) 


where Cyy = B'[yy^] Cyy = E [yy^]- The first inequality follows from the fact that 
{aby < 2{a^ by and the second inequality is shown in Appendix |Bl The first inequality 
splits the MSB into two terms, one contributed by quantization error and the other by channel 
noise. The bound is tight when a is approximately equal to b, that is, when the two MSB 
contributions are approximately the same. Burthermore, the upper bound in (l32l) follows similar 
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to ([7]) and (fTTI) and, thus, are tight under similar eonditions. Similar to the AF ease, even though 
the upper bounds may not be tight in general, they exhibit similar behaviors as the aetual MSE 
with respeet to the sensor densities and transmission thresholds and, thus, are used as objeetive 
funetions in our problem. 

To further analyze ^obj,DF(A, 7) in (l32l) . we note that the elements in Cyy and Cyy ean be 
derived as 


{^yy}hj ~ 


alp{vi,Vj)aiaj, 


for i = j 
for i ^ j, 


(34) 


and 


where 


{^yy} * J 


for i=j 

h^i+^iAqAQtjiqjPijaiaj, for i^j, 


(35) 


Qi = Pr(6i,fc = = 1, Oi = 1), 


(36) 


with 



Efii ^ for i = J 

Ef=i ESi for t ^ j, 


(37) 


T^{i,k},{j,i} = l-2[Pr(6i,fe = 1, bj^i = 0\oi = 1, Oj = 1) + Pr(6i,fc = 0, bj^i = 1 |oi = 1, Oj = 1)] (38) 

for fc = 1,..., i?i, / = 1,..., Bj, i,j = 1,..., M. The derivation of Cyy and a more explieit 
expression of TT{i^k},{j,i} can be found in Appendiees [C] and iDj respeetively. By taking df(A, 7) 
as the objeetive funetion, the seareh for the optimal sensor densities and energy thresholds ean 
be formulated as 


min 

^,7 

^obj,DF(A, 7) 

(39a) 


M 


subjeet to 

5 ^A,A<A, 

i=l 

(39b) 


0 < AjA< eA, for i = 1 ,..., M. 

(39e) 


For the reasons that will be evident later, we perform the ehange of variables where A* is 
replaeed with ai = Fe^{'yi)\iA, for alH = 1,..., M. Then, by omitting the terms that are not 
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relevant to the optimization, the problem ean be written as 


min Jdf(q:, 7) 

(40a) 

a,7 


M 


\ "V Oil -r 

subject to < A, 

(40b) 

0 < a* < exPefyi), for z = 1 ,..., M, 

(40o) 


where 


7) — 


(tr($D„D„))' 


-1 


(41) 


tr($DQ,(Cyy + Cyy)Dc 

Notiee that, in Jdf(ck, 7), the dependenee on a. lies in Dq,, Cyy, and Cyy whereas the 
dependenee on 7 lies only in Cyy. The optimal solution of this problem is still diffieult to find 
due to the non-eonvexity of the problem. However, an approximate solution ean be found by 
using an alternating optimization algorithm ll45l . where a. and 7 are optimized in turn while 
keeping the other fixed and the proeess is repeated iteratively until there is no appreeiable 
deerease in the objeetive funetion. The algorithm is guaranteed to eonverge sinee the objeetive 
is bounded below and is minimized in eaeh step of the algorithm, but may eonverge to only a 
loeal minimum in general. Details of the optimization of a and 7 are deseribed in the following 
subseetion using the Bernoulli energy arrival model as an example. 


A. Solution for the Bernoulli Energy Arrival Case 

Reeall that, in the Bernoulli energy arrival ease, the energy threshold 7* ean be set, without loss 
of generality, as a multiple of 5. In this ease, we have qi = Q (^\J ^^d Fe^{yi) = 

By relaxing the integer eonstraint on 7* into the linear eonstraint 7* > 5, the optimization problem 
beeomes 


min 

q:,7 

7dfK7)- 

(42a) 


M 


subjeet to 

U Vii 

(42b) 


0 < a* < esPibjyi, for z = 1 ,..., M, 

(42e) 


•ji> S, for z = 1 ,..., M, 

(42d) 
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where the term in the numerator is 

M M 

*=1 j¥=i 


2=1 


with 


C^Al^)= .^.. ..M \ 


hhi 


1-4 


-Bi 


and 




Cm (7.7.) = .... .. Ql 


+ Q 




(44) 


hH 


Q 


{2Bi- l){2^j-l) 

for i A h ^rid the term in the denominator is 

4f(«,7) = 5^5^0vC>jja?a|- 


<B, 


Vi,j ) 


(45) 


M M 


(46) 


i=l j=l 

Let and 7^^^ be the solutions obtained in the £-th iteration of the proposed algorithm. In 
the {i + l)-th iteration, we perform the optimization in the following two steps. 

Step 1 (Optimization of ck): In Step I, we first find the optimal value of a given 7 = 


That is, we find 


where 


ct 


(^+1) = argmin Jdf(ck, 7*-^^), 


(47) 


OL^Ti 


(^+1) 




AA - ps 

/ —< A and 0 < Oj < Vz 

tr p‘^ 1? 


(48) 


is the constraint set on o: in iteration £ + 1. By 
be written explicitly as 


, the optimization problem in (l47l) can 


E"iC’"a? + E"iE,,.iA 


mm 

CK 




E m \~^m a 2 
2=1 X/7 = l 


2^,2 
j 


M (^) 

-hrClj < A, 

.=1 


0 < a* < eKPi5/Ai \ for i = 1 ,..., M, 


where - A,iAl + + CiAlT)] and Aj = + 

AAlf\lf'^y\- The optimization problem is nonconvex but can be solved approximately using 




(49a) 


(49b) 


(49c) 
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the condensation method, similar to that in Algorithm [TJ with monomial approximations of the 
denominator of (I49al) . 

Step 2 (Optimization of 7 ): In Step 2, we find the optimal value of 7 given a. = 


i.e., we find 


where 


^(^+1) ^ argmin 7), 


(^+1) 


(50) 




.1 




M (£+1) . 

7 i<A and ( 5 < 7 i<-^,V* 


E 


(51) 


i=l Oil 

is the constraint set on 7 in iteration £+1. Since Jdf(q:, 7 ) depends on 7 only through {Cij}iy=i> 
the optimization over 7 can be formulated equivalently as 


M 


miiii *0,1(7*) l« 


7fc-*^-l' j=l 


,7+i)y 


M 


(«rT 


') 7 • (52) 

i=l j^i 

Notice that the objective function is not a convex function of 7 . However, by applying the upper 
bound Q{u) < 2exp ( —for u > 0, and by taking the high SNR approximation, we have 


^- -e 

Sl,l\ U) rvj 3(2^i ~ 1) 


(53) 


and 


Cijiliy Ij) 


< 


W'^2Bi+Bj 


771,771 

Sx ^ Ba 


{2Bi - - 1) 


Vi,je 


(54) 


for i j. Then, the optimization problem in (1521) can then be approximated as the problem 
below: 


mm 

7 


M 

E 

2 = 1 


0*,*(«f+'))'2(2^» + l) _ A? 
V / ^ ^ 


3(2^* - 1) 


7i 




- ia^ob) ia 


*=1 jV* 


(2^) _ i){2Bj - 1) 


M (£+1) 

subject to ^ ^ ^ 7* < A 


< 7* < for i = 1 ,..., M, 


(55a) 

(55b) 

(55c) 
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Algorithm 2 Alternating Optimization Algorithm of a and 7 

Initialization: Set £ = 0, e > 0 and find a feasible initial solution 

Iteration £ + 1: 

(i) Find by solving (|4^ in Step 1 using the condensation method. 

(ii) Find 7^^+^) by solving (1551) in Step 2. 

(iii) Repeat (i)-(ii) until 

|JDF(a''+‘'.7<'+'>) - /DF(a<'>,7">)l ^ 

JDF(a<'>,7l'>) 

Take and 7^^+^) as the desired solution. 


which is convex and can be solved efficiently using standard numerical approaches, such as the 
interior point method ll42ll . 

By alternating between the optimization problems in Steps 1 and 2 until convergence, the 
desired approximate solution of a and 7 can be obtained. The alternating optimization algorithm 
is summarized in Algorithm |2l 

B. Extension to DF Systems with Parity Check Bits 

In this section, we have so far investigated DF systems where raw bits are transmitted for field 
reconstruction. In practice, the addition of parity bits to the transmitted signal is often considered 
to allow for error detection at the data-gathering node. In this case, the message transmitted by 
the sensor can be treated as an erasure if an error has been detected and as error-free, otherwise. 
The sensor deployment strategy can then be derived similarly in this case. 

Specifically, let us consider the simple case where only an even parity is used. In this case, 
the bit sequence transmitted by a sensor in subregion V* can be written as bj = [bj 6* s.+i], 
where is chosen such that i® ^^en and is referred to as the even parity bit. By 

assuming that an error is always detected and treated as an erasure when it occurs, the effective 
received signal at the fusion center can be written as 


Vi 


( 56 ) 
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Oi = 


(57) 


where 

0, an error is deteeted 
Oj, otherwise. 

Following similar proeedures as in the previous seheme, an MSE upper bound of DF systems 
with a one-bit parity ean be derived as 

^ <crl-tT (^»DaDqCyyDqDa) 

(tr($D^D|))^ 


< 


tr (^DdDqCyyDqD 
where Dq = diag((l — , (1 — and 


— ^obj,PB(A, 7) 


(58) 

(59) 




for i=j, 


a: 


:p{vi,Yj){l-qi)^^+\l-qj)^^+^aiaj, for i^j. 


(60) 


For suffleiently small bit error probability, i.e., for Ri 0 Vi, the objeetive funetion ean be 
approximated as 

{obj,PB ~ -(Ol) 


where 


and 


4b(A,7) 

M 


(62) 


, 1=1 


M M 

Jpb{^, 7 ) = 5^ + (^l + (^IWi + (63) 

1=1 i=l j^i 

By further adopting the Bernoulli energy arrival model, the optimization problem ean be written 
explieitly as 


M M 


max —— 

A,7 ^ 


a. 

^2 2 

i=i j=i '* h 


AS ^ A 2 a 2 

Ec.|+ 1 : 1 : 72 ^ 


(64a) 


i=l 

M 


*=1 j¥=i 


subjeet to ^ A* < A, 

i=l 


(64b) 

0 < Aj < eA < 1 , i = 1 ,. 

..,M, 

(64c) 

7i = Zi5, Zj G M, i = 1,. 

..,M, 

(64d) 
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where Ci = 0*,*(cr^+cr^+crfjp^5^ and Gij = 0ijcr2p(vi, Vj)pjp^j6^. Notice 

that the problem is similar to that obtained in the AF case and, thus, can be solved following 
similar procedures as in Section ITTT-Ai 


V. Simulations and Performance Comparisons 
Appendix A 

Equivalence of the LMMSE Estimator in Q 

Eet Oi = {i\oi = 1} and Oq = {z|oi = 0} be the index sets of the subregions with and 
without a transmitting sensor, respectively, and let |C>i| = r and |C>o| = M — r. Moreover, let 
z = [ 2 : 1 ,..., 2;^]^ be the r X 1 vector obtained by removing the entries in y that correspond to 
the indices in Oq. In this case, = E'[x(v)z^|o] is equivalent to the vector with 

the entries in Oq removed, and 0°^ — -E'[zz^|o] is equivalent to the matrix C°y with the rows 
and columns in Oq removed. The eigenvalue decomposition of 0°^ can be written as 


r 

CL = ^kUkul, (65) 

k=l 

where Pk, for fc = 1,..., r, is the k-th eigenvalue, labelled such that pi > g 2 > ■ ■ ■ > Pr > 0, 
and Ufc = [uk,i, ■ ■ ■ is the corresponding eigenvector. Then, by constructing the vectors 

Vfc = [Vk,i, ..., Vk,MV, for fc = 1,..., r, such that 




0, if Oi = 0, 


( 66 ) 


(i.e., by padding zeros into the vector at locations corresponding to the indices in Oq), we 
have 


r 

Cyy = (67) 

k=l 

Notice that vi,..., are linearly independent and, thus, are eigenvectors corresponding to the 
non-zero eigenvalues pi,..., pr of Cyy. 

Therefore, the EMMSE estimator obtained by using only the received signals r*, for i E Oi, 
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can be written as 


po po-1 

'^a:(v)z'^zz ^ 


(68) 

r r 

i=l j=l 

(69) 

j=i j=i 

f r \ 

'^k,iUk,j 1 Zj 

Kk=l J 

(70) 

M M 

EEfcsw,}*! 

i=l j=l 

^ i?fc 'Pfc.i'Pfc Uj 

(71) 

M M 

EE{c;(v)y}.< 

i=l j=l 

> 

1 

(72) 

po po t,. 


(73) 


where the equality in (TtT]) follows from the fact that Vk,i = 0, Vfc for i G {i|oj = 0}. 

Appendix B 

Derivation of the Inequality in (I^ 

To derive the inequality in (l32l) . let us define the term inside the integral of (l3TI) as 


£(k,v) ^E[|a;(v) -Fy|2]+E[|k'^y-Fy|2] (74) 

= ^X~ 2Ca;(v)yk + k^ (Cyy + Cyy) k (75) 

where C 2 ;(v)y = E[x{v)y'^]. To minimize £(k, v), we set = 0, which yields 

k*(v) = argmin£(k, v) = (Cyy + Cyy)"^C^(v)y^, (76) 

keM*^ 

and the corresponding minimum value is 

£(k*(v), V) = - C.(v)y(Cyy A Cyy)-'C.(v)/ (77) 


where 


Ca;(v)y Cx(v)mDa ~ C3;(v)xDa (78) 

and Cj;(v)m — E[x{\')m^] with m = being the vector of quantized sensor 

measurements. The approximation in (17^ is made by assuming that the quantization error e* is 
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uniformly distributed and is uncorrelated with its input Xi. By substituting (1771) into (|3TI) . the 
upper bound of the average MSB becomes 

2 


Supper,DF ^ |y| j C 3 ;(v)y(Cyy + Cyy) C3;(v)y dv 




approx, DF 5 


where 


0 = 


|V| .Uv 


^x(v)m^a:(v)m<^V 


(79) 

(80) 


(81) 

(82) 


~ l-i;! [ ^x{v)yc^oc{v)yid'V — 

I 7vev 

and $ is defined in (fT3l) . Moreover, using Lemma[T]and the procedure similar to that of (fT5]) - (fT^ . 
the average MSB can be further upper-bounded by 

(tr($D„D,))2 


■Cupper, DF ^ 2 


o-^ - 


tr ($DQ,(Cyy -f Cyy)Dc 


(83) 


Appendix C 

Derivation of {Cyy}i,j: 

Bet us define hi^k — bi^k — bi,k and Oj = {oj = 1} as the event that Oj = 1 for i = 1,, M. 
The diagonal element of Cyy can be written as 

B, 


Em = E 


i.bi,k - bi,k)‘^^‘ 


,fc=l 


Bi 




,fc=i 

Bi 


E 


I2 21 

bi.koAoi 


l-2fc 


Bi 

k=l l^k 
Bi 


E 


bi^kbi^lOj^ |Oj 




(84) 

(85) 


h^khiAni Pr(a)2-'=2-' . (86) 


= 2“-A|, y] E Pr(n,)2-“ + E E ^ 

Vfc=l k=l l^k 

By adopting BPSK modulation with transmit power ei/Bi per bit, the bit error probability of 
the sensor located in V* is 


Qi = Pr(6i^fc = l\bi^k = 0, Oj) = Pr(6i,fc = 0\bi^k = 1, (2*), 


( 87 ) 
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Kk 


n, 


for k = 1,Bi, i = 1,, M. The expectation inside the first term of (18^ is 

= 1 ■ Pr(6i,fc 7^ bi^kl^i) + 0 ■ Pr(6i^fc = Kkl^i) 

= Pribi^k = 0|6i,fc = l,fii)Pr(6i,fc = ll^i) + Pr(&i,fc = = 0, fi*) Pr(6i,fc = 0|fii) 

= Qi- 


The expectation inside the second term of (1861) is 

bi^kbi.^1 1 PT(bi^k 1; bi^i 11 PT(bi^k 1; b-i^i 11 1 ^ 2 ) 

-Pr(6i,fc = 1, \i = -l\ni)-Pr(bi^k = -1, hi = l|^i) 

= Qi [Pribi^k = 1, bi^i = 11 fii) + ^hhk = 0, bi^i = 01fii) 

- Pr(6j,fc = 1, b^ = 0\Pli)-PT{bi^k = 0, b^ = l|fii)] 

= qi[l-2{PT{bi^k = f,bi,i = 0\fli) + PT{bi^k = 0,bi,i = l\ni))], 

where the last equality follows from properties of the natural binary code. By substituting the 
above into (f8^ . we get 


Bi 


EiVi] = '^aiqa + g-r/ija* 


,fc=i 


= n, 


Qi 


1 - 4--®* 


+ qiVi,j 


( 88 ) 

(89) 


The off-diagonal elements in Cyy can be obtained similarly. 


Appendix D 

Evaluation of 7i{i,k},{j,i}- 

To evaluate 7r{i,fc},{j,z} in (l38l) . the joint probabilities Pr(6j^fc = 1, bj^i = 0|f2i, Plj) and Pr(6i^fc = 
Oybj^i = I|f2j,f2j) are required. To do this, let us define 

Ilk ={xi e R\bi^k = z} 


U 


^1,1 —1 




Bi—k — l 


2"'-^ IAq, 


. i=l 


. z=l 


(90) 

(91) 
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as the set of measurement values Xi at sensor i that yields bi^k = 2 ;, where z G {0,1}. The 
approximation follows the faet that Pr(|a:j| > W) ~ 0. Notiee that ean be written as the 
union of disjoint intervals. Therefore, for i = j and k < I, we have 


Pr(6i,fe = = Pr(xi G (1 


(92) 




^i,l >• • • — 1 ) 

^ i.k —-^1 —■^2 


Q 


ELi 


+ 




-Q 


+2^-^ - 2^-1 Aq^ 


(93) 


The result is similar for k > 1. Moreover, sinee Xi and Xj, for i ^ j, are jointly Gaussian 

random variables with mean 0, varianee af = + a^, and eorrelation eoeffieient p(vj, Vj) = 

2 

p(vj, Vj), we have, for i ^ j, 

Pr(6i,fc = l,bj^i = 0|ni,nj) 


= Pr(xi G X/_fc, Xj G I^jj) 


3 


'xj&il Jxi&xl. 2'Kal^Jl- p^{wi,Wj) 


dxidxj 


(a) 


(b) 




TTCri 


f ^ " 


dx'. 


s? 
_ 3 

„ 2cr- 
e a; 




E 


'(EUXj,t2®^ ‘+ 2 -^^ ^)Aq, 


E 


_ ■? 
/n 2cr - 

e 


6 { 0 ,i}A,i =0 e{ 0 ,i}A,fe=i 


Q 


- 2 ^“-^ j Aq,-p(v„v,)x 


Q 


E"=i + 2^'-i Aq- p(v,, v,)5; 


Vf^i(l-P^(vi,Vj)) 
where (a) follows from the ehange of variable 

Xi-p(Vi, Vj)Xj 


dxi 


X, = 


V^l(l-P^(vi,v^)) 


(94) 

(95) 

(96) 


(97) 
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withT^; ^ I 

The probability Pr(6j ^ = 0,6 


G M|6j^fc = 1 } and (b) follows from the definition of and 
j^i = ean be evaluated similarly. 
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